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1  Introduction 

The  South  China  Sea  (SCS),  Yellow/East  China  Sea  (YES),  and  Japan/East  Sea 
(JES)  are  major  east  Asian  marginal  seas  (EAMS).  The  complex  topography 
includes  the  broad  shallows  of  the  Sunda  Shelf  in  the  south/southwest  of  SCS; 
the  continental  shelf  of  the  Asian  landmass  in  the  north,  extending  from  the 
Gulf  of  Tonkin  to  the  YES;  a  deep,  elliptical  shaped  SCS  and  JES  basins,  and 
numerous  reef  islands  and  underwater  plateaus  scattered  throughout  (Fig.  la). 
The  shelf  that  extends  from  the  Gulf  of  Tonkin  to  the  YES  is  consistently  near 
70  m  deep,  and  averages  150  km  in  width. 

The  EAMS  is  subjected  to  a  seasonal  monsoon  system.  From  April  to  Au¬ 
gust,  the  weaker  southwesterly  summer  monsoon  winds  result  in  a  wind  stress  of 
just  over  0.1  N/m2.  From  November  to  March,  the  stronger  northeasterly  winter 
monsoon  winds  corresponds  to  a  maximum  wind  stress  of  nearly  0.3  N/m2.  Re¬ 
cent  observational  studies  show  that  the  EAMS  is  energetic  and  has  multi-eddy 
structure.  For  example,  the  SCS  synoptic  eddy  structure  was  identified  in  May 
1995  using  the  airborne  expendable  bathythermograph  (AXBT)  data  (Chu  et 
al.,  1998a),  the  eddy  spatiotemporal  scales  in  the  YES  were  identified  using  the 
Navy’s  Master  Oceanographic  Observational  Data  Set  (MOODS)  during  1929- 
1991  (Chu  et  al.,  1997a, b),  and  the  seasonal  JES  multi-eddy  structure  from  a 
composite  analysis  on  the  U.S.  National  Centers  for  Environmental  Prediction 
(NCEP)  monthly  SST  fields  during  1981-1994  (Chu  et  al.  1998b). 

The  fundamental  scientific  issues  axe:  What  axe  the  dynamical  balances 
contxolling  mesoscale  eddy  variations  in  the  EAMS,  including  the  effect  of  wind 
and  thermohaline  forcing,  wave  processes,  current  instabilities,  coastline  geom¬ 
etry  and  topographic  slope?  What  is  the  role  of  coastal  air-ocean  coupling 
on  oceanic  and  atmospheric  dynamics  of  mesoscale  eddies?  To  study  these 
problems,  a  coastal  atmosphere-ocean  coupled  system  (CAOCS)  was  developed 
at  the  Naval  Postgraduate  School.  The  model  domain  (67°-142°E,  7°S-55°N) 
covers  the  whole  EAMS  and  surrounding  land  and  islands.  The  surface  fluxes 
of  water,  heat  (excluding  solar  radiation),  and  momentum  axe  applied  syn- 
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chronously  with  opposite  signs  in  the  atmosphere  and  ocean.  Flux  adjustments 
are  not  used. 

2  CAOCS  for  the  EAMS  Nowcast /Forecast 

2.1  Ocean  Component 

2.1.1  Numerics 

The  CAOCS  ocean  component  (i.e.,  the  POM)  is  a  time  dependent,  primitive 
equation  circulation  model  on  a  three  dimensional  grid  that  includes  realistic  to¬ 
pography  and  a  free  surface.  From  a  series  of  numerical  experiments  (Chu  et  al., 
1998a;  1999a, b;  2000b),  the  qualitative  and  quantitative  effects  of  non-linearity, 
wind  forcing  and  lateral  boundary  transport  on  the  EAMS  are  analyzed,  yielding 
considerable  insight  into  the  external  factors  affecting  the  regional  oceanogra¬ 
phy.  We  use  a  rectilinear  grid  with  horizontal  spacing  of  0.25°  by  0.25°  and  23 
nonuniform  vertical  a  levels.  The  model  uses  realistic  bathymetry  data  from 
the  Naval  Oceanographic  Office  Digital  Bathymetry  Data  Base  with  5  minute 
resolution  (DBDB5). 

2.1.2  Open  Boundary  Conditions 

Closed  lateral  boundaries,  i.e.,  the  modeled  ocean  bordered  by  land,  were  de¬ 
fined  using  a  free  slip  condition  for  velocity  and  a  zero  gradient  condition  for 
temperature  and  salinity.  No  advective  or  diffusive  heat,  salt  or  velocity  fluxes 
occur  through  these  boundaries. 

At  open  boundaries,  we  use  monthly  varying  volume  transports  from  a  global 
inverse  model  (Chu  and  Fan,  2000)  for  the  open  boundaries  (Fig.  la).  When 
the  water  flows  into  the  model  domain,  temperature  and  salinity  at  the  open 
boundary  are  likewise  prescribed  from  the  climatological  data.  When  water 
flows  out  of  the  domain,  the  radiation  condition  was  applied, 

^(0,S)  +  I/«^(0,S)  =  O  (1) 

where  the  subscript  n  is  the  direction  normal  to  the  boundary. 

2.1.3  Initial  Conditions  and  Model  Initialization 

Before  coupling  to  the  atmospheric  model  that  is  MM5,  the  POM  model  was 
integrated  for  four  years  and  four  months  from  zero  velocity,  and  January  cli¬ 
matological  temperature  and  salinity  fields  and  forced  by  monthly  mean  surface 
wind  stress  from  the  Comprehensive  Ocean  and  Atmosphere  Data  Set  (COADS) 
and  by  the  restoring-type  surface  heat  and  salinity  fluxes  which  are  relaxed  to 
the  surface  monthly  values.  The  final  states  of  temperature,  salinity,  and  veloc¬ 
ity  were  taken  as  the  initial  ocean  conditions  for  1  April  1998. 
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2.2  Atmospheric  Component 

2.2.1  Model  Description 

We  take  Version-3  of  the  Fifth-Generation  Penn  State  University/NCAR  Meso- 
cale  Model  (MM5)  as  the  atmospheric  component  for  CAOCS.  The  MM5  is  a 
limited-area,  nonhydrostatic,  terrain-following  sigma-coordinate  model  designed 
to  simulate  or  predict  mesoscale  and  regional-scale  atmospheric  circulation.  It 
has  been  developed  at  Penn  State  and  NCAR  as  a  community  mesoscale  model 
and  is  continuously  being  improved  by  contributions  from  users  at  several  uni¬ 
versities  and  government  laboratories. 

2.2.2  Land  Surface  Parameterization 

The  horizontal  resolution  is  60  km.  Sixteen  pressure  levels  were  used  with  the 
top  at  10  mb.  A  soil  water  availability  function  was  used  for  computing  the  soil 
water  content  by  using  10  specified  vegetation  types.  Dominant  deciduous  and 
coniferous  forests  cover  the  most  land  areas.  There  are  small  areas  of  tropical 
forest  in  Taiwan,  and  Indo-China  Peninsula. 

2.2.3  Initial  and  Lateral  Boundary  Conditions 

The  initial  and  horizontal  lateral  boundary  conditions  for  wind,  temperature, 
water  vapor,  and  surface  pressure  are  interpolated  from  analyses  of  observations 
from  European  Center  for  Medium-Range  Weather  Forecasts  (ECMWF).  For 
the  present  study,  they  were  projected  on  a  spectral  T42  grid,  and  a  vertical 
resolution  of  14  pressure  levels.  The  initial  conditions  were  the  fields  on  00Z 
April  1,  1998.  The  lateral  boundary  conditions  were  provided  via  a  relaxation 
method  and  updated  every  12  hours. 

2.3  Ocean-Atmosphere  Coupling 

Atmospheric  and  oceanic  surface  fluxes  of  water,  heat  (excluding  solar  radia¬ 
tion),  and  momentum  are  of  opposite  sign  and  axe  applied  synchronously.  Flux 
adjustments  are  not  used.  We  use  the  same  flux  parameterization  as  the  MM5 
for  CAOCS.  The  only  difference  between  the  two  is  the  use  of  sea  surface  tem¬ 
perature  (SST)  data.  In  stand-alone  MM5  model,  SST  is  prescribed  as  a  given 
parameter.  However,  in  the  CAOCS,  SST  is  predicted  by  the  POM  model  at 
each  time  step  for  the  atmospheric  component  (MM5). 

3  Predicted  EAMS  Circulation 

We  integrate  the  CAOCS  model  for  four  months  from  1  April  1998.  Figures 
lb,c,d  show  the  predicted  temperature  and  velocity  vector  fields  at  three  dif¬ 
ferent  depths  (surface,  50m,  and  200  m)  in  June  1998.  The  North  Equatorial 
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Current  (NEC)  was  well  predicted  in  all  depths.  After  encountering  the  west¬ 
ern  boundary  along  the  Philippine  coast,  the  NEC  bifurcates  into  the  northward 
flowing  Kuroshio  and  the  southward  flowing  Mindanao  Current  (MC). 

The  MC  and  the  simulated  New  Guinea  Coastal  Undercurrent  flow  equa- 
torward,  feeding  both  the  North  Equatorial  Counter  Current  (NECC)  and  the 
flow  from  the  Pacific  to  the  Indian  Ocean,  that  is,  the  Indonesian  throughflow. 
The  MC  is  an  essential  part  of  the  shallow,  wind  and  buoyancy  flux  driven  sub¬ 
tropical  circulation  cell  that  carries  the  subtropical  waters  to  the  equator.  It 
is  typically  colder  and  fresher  than  the  South  Equatorial  Current  (SEC)  and 
the  New  Guinea  Coastal  Current  over  the  thermocline,  due  to  ventilation  in  the 
North  Pacific  subtropical  gyre.  The  northward  flowing  Kuroshio  is  the  major 
western  boundary  current  in  the  Pacific  Ocean.  The  Kuroshio  Water  intrudes 
into  the  SCS  through  the  southern  Luzon  Strait  (the  strait  between  Taiwn  and 
Phillippines),  loops  around  anticyclonically,  bifurcates  southwest  of  Taiwan  into 
the  northeast  flowing  South  China  Sea  Warm  Current  (SCSWC)  and  southeast 
flowing  current  which  exits  the  SCS  through  the  northern  Luzon  Strait.  In  the 
SCS,  a  multi-eddy  structure  was  predicted  with  central  SCS  warm  anticylonic 
eddy  with  several  surrounding  cyclonic  eddies.  The  continuation  of  the  SCSWC 
becomes  the  Taiwan  Warm  Current  after  passing  the  Taiwan  Strait.  Major  cur¬ 
rent  system  in  the  EAMS  were  realistically  predicted  using  the  CAOCS.  The 
flow  pattern  is  consistent  with  recent  observational  studies  (e.g.,  Chu  and  Li 
2000)  and  earlier  studies.  Interested  readers  are  referred  to  a  recent  review 
paper  on  the  EAMS  (Su  1998). 

The  CAOCS  model  was  assessed  using  the  data  collected  from  the  inten¬ 
sive  observational  period  (IOP)  of  the  international  South  China  Sea  Monsoon 
Experiment  (SCSMEX)  and  from  the  Navy’s  AXBT  measurements  in  April- 
August  1998. 

4  CONCLUSIONS 

The  CAOCS  has  capability  to  predict  the  current  system  and  thermal  struc¬ 
ture  of  the  EAMS.  Comparing  to  the  SCSMEX  and  AXBT  measurements,  the 
CAOCS  successfully  predicted  the  South  China  Sea  multi-eddy  structure  during 
the  SCSMEX  IOP  in  1998.  These  warm-core  and  cool-core  eddies  have  radii 
varying  from  100  km  to  300  km  and  maximum  tangential  velocities  ranging  from 
10  cm/s  to  20  cm/s.  The  cool-core  eddies  are  cyclonic  and  the  warm-core  eddies 
are  anticy clonic. 
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Fig.  1.  Predicted  temperature  and  velocity  fields  in  June  1998 
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